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We study the BCS-Bose Einstein Condensation (BEC) crossover of a three dimensional spin 
polarized Fermi gas with Rashba spin-orbital-coupling (SOC). At finite temperature, the effects of 
non-condensed pairs due to the thermal excitation are considered based on the GqG pair fluctuation 
theory. These fluctuations generate a pseudogap even persistent above Tc- Within this framework, 
the Sarma state or the spin polarized superfiuid state and polarized pseudogap state are explored 
in detail. The resulting Tc curves show that the enhancement of pairing due to the SOC roughly 
cancels out the suppression of pairing due to the population imbalance. Thus we observed that in 
a large portion of the parameter space, the polarized superfiuid state are stabilized by the SOC. 


I. INTRODUCTION 

Ultracold Fermi gases with tunable interactions has 
been the focus of a lot of experimental and theoretical 
works. The advantage of cold Fermi gases is that they are 
much more controllable than other ordinary condensed 
matter materials. The interaction between the fermions 
usually has very short effective range and can be char¬ 
acterized by one parameter, the s-wave scattering length 
a, or more conveniently, dimensionless scattering length 
l/{kpa), with the Fermi momentum kp. By applying a 
external magnetic field, one can tune the Zeeman split¬ 
ting energy between the bound state of closed channel 
and the continuum threshold of the open chann^ which 
will lead to the so-called Fashbach resonance [J-Q when 
the above two energy levels line up. This phenomenon 
allows one to tune 1/ {kpa) from very negative to a large 
positive number. Correspondingly, the interatomic in¬ 
teraction varies from a weak attraction to a very strong 
attraction. Therefore, the cold Fermi gases provide an 
experimental realization of the early theoretical ideas of 
BCS-BEC crossover[l-0 and are also very useful for test¬ 
ing the many-body theories. 

A natural interesting point in the crossover is the 
so-called unitary limit when the scattering length di¬ 
verges. The unitary limit of Fermi gases is intrinsi¬ 
cally strong correlated and there is no small parame¬ 
ter to expand with. The BCS-BEC crossover theory de¬ 
scribes the loosely bound Cooper pairs evolve into tightly 
bound bosonic pairs with increasing attractive interac¬ 
tion strength. It captures the two weakly interacting lim¬ 
its and can also give quantitative account of the unitary 
limit in between. When considering the finite tempera¬ 
ture physics of the crossover, the contribution of the non- 
condensed pairs become important due to the stronger 
than BCS attraction and these effects must be consid¬ 
ered in a way consistent with the BCS ground state. 

Soon after the experimental evidence of the fermionic 
superfluidity was achieved, one was able to adjust the 
number density of spin up and spin down particles. Since 
the singlet pairing requies equal number of spin up and 
down particles, this spin polarized Fermi gases generated 


a lot of possible novel phases Q, such as Sarma or polar¬ 
ized superfiuid phase p. polarized pseudogap phase, 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase [ill- 13. 
or the phase separations between different phases [l4[l^. 
The stability of these phases and the transitions between 
them have also been studied by many authors [iH. IT^. 

In recent years, another important breakthrough in ex¬ 
periments is the realization of the synthetic non-Abelian 
gauge field and the spin-orbital-coupling (SOC) in cold 
atomic gases [l^. By applying two counter propagating 
Raman laser beams and a transverse Zeeman field to Bo¬ 
son atoms with multiple components, a synthetic SU(2) 
non-Abelian gauge field can be generated by appropriate 
choice of the laser frequency and Zeeman splitting. One 
can show that the low energy effective theory contains 
a Raman type of SOC terms. Similar scheme can also 
be applied to the fermion atoms, and the SOC of Fermi 
gases has already been generated in experiments (T^. . 

The SOC of Fermi gases opens up the possibilities to ex¬ 
plore many novel physics in cold atoms. It also provides 
us yet another way to achieve the BCS-BEC crossover. 
As pointed out in [2l[, with Rashba type of SOC, there 
appears a new type of two-body bound state of fermion 
pairs even for a < 0, while two-body bound state does 
not exit in the BCS side without SOC. This suggests 
that the pairing strength between fermion pairs are en¬ 
hanced by the Rashba SOC term. Therefore, increasing 
the SOC coupling is equivalent to push the Fermi gases to 
its deep BEC limit, thus is another way to get BCS-BEC 
crossover. 

There already appears a lot of theoretical works 
on the SOC Fermi gases with or without population 
imbalance [^. However, other than a few exceptions!^, 
most of these works are based on the mean field theory 
which should be more appropriate for the deep BCS or 
very low temperature cases. The reason is that, as we 
mentioned before, in the unitary limit or BEC side of 
the crossover, the attraction is strong enough to support 
two-body bound pairs. At the finite temperature, there 
will be substantial amount of non-condensed pairs which 
generate an energy gap in single particle spectrum. In 
contrast to the superfiuid order parameter, this pseudo¬ 
gap does not signal symmetry breaking. Thus it can be 
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persistent above which leads to a non-Fermi liquid 
normal state. Different pair fluctuation theories differ¬ 
entiate themselves from each other in the detailed form 
of pair propagator or T-matrix (23 - l^ . In this paper, 
we follow the so-called GqG pair fluctuation theory 
which is pa rtly inspired by the early work of Kadanoff 
and Martin[^. The advantage of this theory is its con¬ 
sistency with the BCS ground state and its numerical cal- 
culability. With both SOC and spin polarization, there 
maybe emerges a lot of exotic phases. To map out the 
whole phase diagram will be a enormous task. In this 
paper, we only consider the simplest Sarma state or po¬ 
larized superfluid state and polarized pseudogap state. 
To avoid other possible phases, we will mostly confine 
ourself in the parameter space where the polarized su¬ 
perfluid state is stable. 

This paper is organized as follows. In section |nl we 
present the mean field theory of population imbalanced 
Fermi gases with Rashba SOC. Then we generalized the 
mean field theory to include the GqG pair fluctuation 
effects in section Ell In section HVl we present the nu¬ 
merical results and discussion bases on the previous the¬ 
oretical formalism. We conclude in section 0 


imaginary-time Green’s function can be written as: 


= ’G(tx) 

F{t,x.) G{t,x.) J ’ 


where Tr is the time order operator. After Fourier trans¬ 
formation, we find the Green’s function in frequency- 
momentum space: 


g{K) 


G{K) F{K)- 
F{K) G{K) J ’ 


( 3 ) 


Here K = (ia;„,k), = (2n -I- l)7rT is the Matsubara 

frequency for fermion. And its matrix elements also sat¬ 
isfy the following relations 


G{K) = -G{-Kf (4) 

F{K) = -Fi-Kf (5) 


II. MEAN FIELD THEORY 


In this section, we consider the BGS mean field the¬ 
ory of population imbalanced Fermi gases with Rashba 
type of SOG. We will see that the gap equation can be 
rewritten in the T-matrix form, according to the Thou- 
less criterion. Through this form, it is easy to generalize 
the BGS mean field theory to include the GqG pairing 
fluctuation effects. 

The system of two-component population imbalanced 
Fermi gases with Rashba SOG across a Feshbash reso¬ 
nance can be described by the Hamiltonian: 

H = J fj,-5fia^+'Hso]i>{x) 

+g J d^xV’|(x)V'j(x)?/'i(x)V't(x), (1) 


where = {ip^{x),ipj^{x)) is the creation operator, 

/i = + At 4.)/2 is the average of the chemical poten¬ 
tials, dfi = — /i 4.)/2 is the chemical potential differ¬ 

ence, Uso = —i^icfxdx + cTydy) IS the Rashba SOC term, 
ax,(7y,(7z the three Pauli matrices, g < 0 is the bare 
interaction strength of s-wave attractive interaction. In 
this paper, we take ft = fcs = I for convenience. 

It is more convenient to introduce the Nambu spinor 
'k'l' = (^|, 'i/jj, tpi). Because the SOC interaction term 
explicitly depends on the spin indices, the Nambu spinor 
has 4 components, mean field Hamiltonian and Green’s 
functions in Nambu space are 4 by 4 matrices. Then the 


Next we will derive the gap equation, number density 
equation, number density difference equation based on 
the mean field approximation Age = g{'tpx{T, x)'!/)'|-(r, x)). 
In this paper, we only consider the pairing in the spin 
singlet channel. For convenience, we take Age to be real, 

i'6. Ajjp = Age. 

By rewriting the Hamiltonian Eq.(IT|) in the Nambu 
space and making the mean field approximation to the 
two-body interaction[2^, we find that the inverse BCS 
propagator in frequency-momentum space can be ex¬ 
pressed as 


g-\K) 


G^\K) -lAscay' 
iAsc<Jy G^\K) \ ’ 


( 6 ) 


where 


Gq ^{K) = iuJn - (^k - Sga^) - X{K(Tx -b kyay), 

Gq ( at ) = iujji X{kx^x kyCiy'j^ 

with = k^/2m — g. 

A simple matrix inversion gives the following full BCS 
Green’s function in Nambu space 


g{K) 


■G{K) F{K)- 
_ F{K) G{K) \ ’ 


( 7 ) 


where G{K) and F{K) are 2 by 2 matrices with the fob 
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lowing matrix elements 

Gii = ^|a7Kk»7k-<^M(AL + Ck)-7<^MCk£^k] 

1 


7 , 0 ! 


Gi 2 

G 21 

G 22 


+Pk£^k + 7 Pk(Ck — -7777^-J= 77 Ti 

J Ap]^E^{iu:n - 
^ ^ iKkx - iky){a^l. + a7Ck£^g + Pk) 


7,0! 


= E 


4/0k£^k(*‘^" - 7-E'k) 
a7Ck£^k 


^\{kx + iky){a^j 


Pk) 


7 , 0 ! 


4/Ok^^k “ 7-E'k) 

E {«7Kk»7k + + Ck) + 7<^MCkK] 

7 , 0 ! 

+Pk£^k + 7Pk(Ck + - 7777 T- EoE’ 

J 4pkE^{zUJn - lE^) 
a'yXikx - iky)A{^ + Sp) 


= E-' 


4pE^{iujn - 'yE° 


"12 


= E 


jA{p + aSp^) — aASpE^ 


4pE^{iu}n - -fE^) 


= E- 


7 A(p + aSp^) + aASpE^ 

4pE^{iuJn - lE^) 


p„. 


= E 


01,7 


ajXjkx - iky)A{^ - 5p) 

4p-E“(iw„ - 7 £^k) 


and the other two tilted Green’s functions are given by 
Eq.(|4|) and ([5]). Here for convenience, we introduce 7?k = 

^X^{kl + fc2) + Sp'^, pk = The qusi- 

particle energy is given by E^ = ^/^+A^~^r^+2a^ 
and a, 7 = ±. 

With the full BCS Green’s function, we find the self- 
consistent gap equation by evaluating the expectation 
(' 04 ,(r, x)'!/)-|-(t, x)) as follows 

^ _ 

2 
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= -jTr^*a,F(if) 


K 




( 8 ) 


Here J2k = TJ2nJ2k and f{x) = l/[exp(a;/r) 1 ] is 
the Fermi-Dirac distribution function. 

Similarly, the number density equation is 


n = n-f + ni = Tr G{K) 


f) 7^^ 

1 + ^[ 1 _ 2 /(K)] 


k,a 


dp 


(9) 


and number density difference equation is given by 
Sn = — ni = TT^^iazG{K) 


K 




2 ^ dSp 

k,a 


( 10 ) 


It should be noted that the above three equations can also 
be obtained directly by taking derivatives of the thermo¬ 
dynamic potential as in[^. 

We have derived the gap equation from standard mean 
field theory calculations. The gap equation can also be 
obtained from a T-matrix point of view. The T-matrix 
can be thought as the proper pair propagator with ex¬ 
ternal legs been amputated. According to the Thouless 
criterion (31|. the divergence of many-body T-matrix sig¬ 
nals the instability of normal ground state, which leads 
to the condensation of fermion pairs and the formation 
of the superfluid ground state. The many-body T-matrix 
can usually be approximated by the summation of a se¬ 
ries of ladder diagrams as 


t{Q) = 


( 11 ) 


1 + 9X{Q) ’ 

where xiQ) is the mean-held pair susceptibility: 


xiQ) = -Tr^GoiK -Q)iayG{K)i(7y 


( 12 ) 


K 


with Q = (*fl„,q), Xln = 2mTT is the Matsubara fre¬ 
quency of bosons. We note that the pair susceptibil¬ 
ity is made by one full and one bare Green’s functions. 
As pointed by Kadanoff and Martin long time ago[2^, 
this is important to recover the BGS gap equation. The 
Thouless criterion require that the many-body T-matrix 
diverges at Q = 0 , i.e. t“^( 0 ) = 0 , or in explicitly form: 

1 + 5X(0)=0, T<Te. (13) 

It can be checked that equation Ea. (fT3ll is exactly the 
same as the gap equation Eq.®. 

In the same time, the mean-held self-energy can also be 
expressed in terms the T-matrix. According to the Dyson 
equation, the self-energy is dehned as the difference of 
Go (AT) and G{K) as follows 

Esc{K) = Go\K)-G-^{K) 

= -^lMvGa{K)i(jy. (14) 

In the superhuid phase, the T-matrix diverges at Q = 
0 , which allow us to approximate it by a delta function 
located at zero momentum. Therefore, we introduce the 
T-matrix for the condensate pairs 

G7Q) = -Al^5{Ql (15) 

Here the delta function of 4-momentum is 5{Q) = 
(l/T)(5ji_o'5*'^H‘l)i then the self-energy can be rewritten 
as: 


Esc{K) = '^tsc{Q)i(JyGo{K - Q)iay. (16) 
Q 

It seems that it is rather cumbersome to write the self¬ 
energy in such a form, but one can see that in this form, 
the full Green’s function can be treated as the bare prop¬ 
agator dressed by the condensate T-matrix. 
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We can see that only the contribution from the con¬ 
densed pairs is considered in the BCS theory. In the 
BCS-BEC crossover, the attractive interaction gets much 
stronger when one approaches unitary limit or the BEC 
side, therefore one must also include the contribution 
from non-condensed pairs or pairing fluctuations. In con¬ 
trast to the condensed pairs whose total momentum is 
zero, the non-condensed pairs have non-zero center of 
mass momentum and its effects become important if we 
want to consider the finite temperature physics. In the 
following section, we will consider the GqG pairing fluc¬ 
tuation theory, which allow us to generalize the above 
mean field theory results to the finite temperature with 
certain strong correlation effects included. 


III. PAIR FLUCTUATION FORMALISM 

In the mean-field theory, only the contribution from 
condensed pairs was taken into account. We also want to 
include the contribution of non-condensed pairs at finite 
temperature. In general, the non-condensed pairs can 
be described by the amputated pair propagator or T- 
matrix with non-zero total momentum. Here we use the 
simple ladder diagrams to approximate the T-matrix. In 
momentum space, the T-matrix of non-condensed pairs 
are given by 


(17) 

In order to be consistent with the BCS gap equation , we 
define the pair susceptibility as 

xiQ) = ^Tr^ Go(if - Q)iayG{K)zay (18) 

as already discussed in [23, Ull • The above T-matrix may 
looks exactly the same as the mean field case discussed in 
the section m but here the full Green’s function G also 
contains the self-energy dressed by the non-condensed 
pairs which is much more complicated than the BCS full 
Green’s function used in the section m 

The total T-matrix combining the condensed and non- 
condensed contribution is given by 

t{Q) = tUQ) + tpg{Q) (19) 

which give rise the following self energy 

T,{K) = t{Q)iayGo{K - Q)iay 

Q 

= S,,(A) + Epg(A), (20) 

where 

E,,(A) = -AlYayGoiK)iay, (21) 

^pgiK) = tpgiQ)iayGoiK - Q)iay, (22) 
Q 


In fact, the Eg. lfTTl) and (l20l) dehne a set of self-consistent 
equation system for the the unknown function T,pg. With 
some initial guess of Epg, one can determine the full 
Green’s function and further compute tpg. The T-matrix 
tpg can be used to determine Epg again, which closed the 
loop. One can keep the above iteration until a conver¬ 
gence is reached. Once Y,pg is obtained, the gap equation 
can be derived from generalized Thouless criterion 

t-y0) = i + x(0) = 0 (23) 

9 

However, it too difficult to solve the full GqG theory nu¬ 
merically. 

In order to simplify the above mentioned numerical 
procedure, we have to make some further approxima¬ 
tions. Note from the Thouless criterion that tpg is highly 
peaked around Q = 0, hence the self-energy can be ap¬ 
proximated as: 

Epg « -AlgmyGo{K)iay, T < Tc, (24) 

where the pseudogap is defined as: 

Apg = — tpg{Q). (25) 

Q/o 


After this approximation, the Epg takes the same form as 
the BGS self-energy, which greatly simplihes calculations. 
Then the total self-energy can also be written in the same 
form as BGS self-energy, 

E{K) = -A‘^iayGo{K)iay, (26) 

with -I- Apg . This also agrees with the intuitive 

picture that the correction of single particle propagator 
comes from both condensed and non-condensed pairs. 

Since the total self-energy takes the same form as the 
BGS thoery, we expect that the gap equation, number 
equation, number difference equation also take the same 
form as in the BGS mean field theory, except that the 
order parameter Age is replaced by total energy gap A: 


A 

9 

n 


5n 


k,Q; 


(27) 

(28) 
(29) 


To find out the pseudogap, we can compute the pair 
susceptibility by substituting the full Green’s function 
G{K) and free Green’s function Go (AT) into Ea. dT^ . 
Note that G{K) also take the same form as in BGS theory 
but with Age replaced by A. The result is 
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x{Q) = 


■S7(«gk + + pk)[A^(fcg + fcg) - X'^{kxqx + kyQy)] 

4pki^k^k-q 

?7k-q(pk£'£ + 7PkCk + a7fk??k) + S(5//2[7Pk + a7(A2 + + a^k-Eg^] 


E 

s,7,a=it k 


i:[ 


4pk^^k7k-c 


/(e^-q)-/(7ig£) 

90 - 7^^k + ^1 


(30) 


k-q 


where = -^k-q + S7k-q- 

Then the pseudogap can be determined by Eg. (ESI) . 
In order to further simplify the calculations, we can take 
the advantage that the T-matrix is highly peaked around 
(5 = 0. Thus we expand the inverse T-matrix t~g(Q) 
around Q = 0 and only keep the leading terms. The 
gap equation Eg. (1231) ensures that constant term of this 
expansion is zero in the superfluid phase. Therefore, we 
find the following expansion 

^pgiQ) ~ ao(if^n-^q) (31) 

3 

2=1 


Here we introduce the pair dispersion Hq. oq is the in¬ 
verse of spectral weight of fermion pairs and Bi is the 
inverse of the effective mass of fermion pairs along direc¬ 
tion i with i = x,y, z. We have: 


flo 


dxiQ) 


dqo 


Q=0 


o-oBi 


1 d'^xiQ) 

2 dqf 


Q=0 


where i = x,y, z. The detailed expressions of oq and Bi 
are very complicated, which are written out explicitly in 
Appendix O 

Inserting the above expansion to Eg. (1251) . we find 


Ka = (32) 


Here b{x) = I/[exp(a:/T) — 1] is the Bose distribution 
function. Solving Eo. (l27)) - (|29)) will determine the gap 
A, chemical potential y, chemical potential difference Sfi 
below Tc- With these results, the pseudogap Apg can 
be determined by Eg. (15^ and the superfluid order pa¬ 
rameter is given by A^^ = A^ — A^g. The vanishing of 
Asc determines the critical temperature T^, which is also 
equivalent to the condition A = Apg at T^- 

Above Tc, there is no symmetry breaking, thus gap 
equation Eo. (l23l) is not valid any more. We can extend 
our theory from below to above by insisting that 
in the normal phase the energy gap is solely generated 
by the pairing fluctuations. Thus we assume that A = 


Apg for T > Tc- Since Eo. d^ does not hold any more, 
the expansion of tpg around Q = 0 will give a non-zero 
constant term, which can be treat as the effective pair 
chemical potential. 


tpniQ) 


ao 


i^n (^q 


(33) 


The pair chemical potential jib is determined by 

1 


aofJ-b = - + x( 0 ) 
1 


4A 


f) fpOi 

E^[l-2/(^k)] (34) 


This equation will serve as the gap equation above 
Due to non-zero the pseudogap equation above Tc is 
slightly modified as 

(35) 

u-O 

q 

The number density and number density difference equa¬ 
tion Eo. (l28l) . (l2^ are the same as before. Combining the 
number equations with the new gap and pseudogap equa¬ 
tions Eo. (l34|) and Eo. (|M)) . one can determine Sy, /if,, 
and Apg above T^. Therefore we have generalized the 
pair fluctuation theory to the polarized Fermi gases with 
SOC. In the next section, we will show the numerical 
results of the above theory and also discuss its physical 
meaning. 


IV. NUMERICAL RESULTS AND 
DISCUSSIONS 


In this section, we present our numerical results of the 
polarized Fermi gases with Rashba SOC. We will mainly 
focus on the polarized superfluid phase and the polarized 
pseudogap phase, and ignore the possibility of more ex¬ 
otic LOFF phase and phase separations. In the numer¬ 
ical calculations, the bare interaction strength g is re¬ 
placed by the experimental measurable scattering length 
a. These two quantities are related by the following reg¬ 
ularization equation. 


1 TTl V ^ 1 

g Ana ^ 2ek' 


(36) 
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(a) \ =Q.2bvp, 5n=0.1, l/fc^a = -0.5 



(b) \ =Q.2bVp, 5n=0.4, l/kpa = -Q.b 


A/ep 

Apg/^F 

^sc/ 


/ N 




(d) X=Vp, 5n=0A, l/kpa = —0.5 



0.0 0.1 
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T/ep 
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0.2 0.3 

T/ep 
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Figure 1: The total gap A, superfluid gap Asc and pseudogap Apg as functions of T/Ep with interaction l/kpa = —0.5 which 
corresponds to the shallow BCS case. In panel (a),(b) and (c),(d), the SOC coupling constant X/vp = 0.25, 1 respectively. In 
panel (a),(c) and (d),(d), the number density difference 5n/n = 0.1, 0.4 respectively. 


For convenience, we take the Fermi energy and momen¬ 
tum to be unity, which corresponds to taking the density 
as n = I/Stt^. We treat l/^kpa) as the effective interac¬ 
tion. 

Due to the strong attractive interaction, the appear¬ 
ance of pseudogap is quite general phenomena in the 
BCS-BEC crossover of Fermi gases. The advantage of our 
GqG pair fluctuation theory is that the pseudogap effect 
is relatively easy to compute. In FiglTJ we first show the 
behaviors of total energy gap A, superfluid gap Asc and 
pseudogap Apg depending on the temperature. All the 
plots in Fig. [T] have similar features as the Fermi gases 
without the SOC and spin polarization. The total gap 
decrease with the increasing T as in the mean field theory. 
Apg is zero at T = 0. Then, the pseudogap increases with 
increasing T roughly as a power law function, reflecting 
the increasing of thermally excited non-condensed pairs. 
Finally, the pseudogap saturates the whole energy gap, 
which determines the transition temperature Tc- Above 
Tc, the appearance of nonzero /if, makes the pair number 
density become smaller, therefore the pseudogap gradu¬ 
ally dies off. The polarized pseudogap phase is a special 


normal phase with single particle energy gap, which is 
clearly in contrast to the normal Fermi liquid. However, 
we should mention that above Tc, the pseudogap decays 
quite slow and will extend to even above Ep. But one 
should expect, there exist a crossover temperature T* 
above which the system should go back to Fermi liquid 
or Fermi gas. Therefore, the high T part of the above 
plots only has qualitative meanings. 

In the unitary limit and BEC side, the pairing effects 
are very strong. Thus the suppression of the superfluid¬ 
ity by the population imbalance is not very prominent, 
we can expect that in these regime the polarized super¬ 
fluid phase is stable. On the other hand, in the deep 
BCS limit, one should expect that the polarized super¬ 
fluid phase will be replaced by phase separation between 
balanced superfluid phase and normal phase of major¬ 
ity species. In FiglU we shows the results for the shal¬ 
low BCS case l/{kpa) = —0.5, which is roughly on the 
boundary between the stable and unstable polarized su¬ 
perfluid phase. Without the SOC, the pseudogap will 
become non-monotonic as a function of T for the same 
set of parameters of Ijkpa and Xmjn. In that case, there 
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Figure 2: (a) Tc as a function of Ijkpa, with fixed SOC coul- 
ing A = 0.25vf and three different number density difference 
Sn/n — 0, 0.1, 0.3. (b) Tc as a function of l/kpa, with fixed 
number density difference Sn/n — 0.1, with and three differ¬ 
ent SOC coupling X/vf = 0, 0.25, 1. 


exists an upper and lower T^, which clear signals that 
the polarized superfluid state are un-physical and there 
must be a mixture of superfluid state and normal state. 
In FiglU we see that the introducing of SOC term in¬ 
creases the pairing effects and make the polarized su¬ 
perfluid state stable. Even when 6n/n = 0.4, the SOC 
coupling A = 0.25ui^’ is still strong enough to stabilise the 
superfluid state. Comparing panel (c) and (d) to panel 

(a) and (b), we see that the increasing of SOC strength 
generally increase the overall magnitude of pairing gaps. 

In Fig 121 we display the Tc as a function of the inter¬ 
action l/kpa, with different number density difference 
in panel (a) and with different SOC strength A in panel 

(b) . For the balanced case with SOC, we reproduced 
the same result as in Comparing with the balanced 
Fermi gases without SOC, one can see that SOC term did 
not make too much qualitative change to the Tc curve. 
The maximum of Tc is reached around unitary limit, and 
there is a small local minimum located at the near BEC 
side. These features are quite general in the GqG pair 
fluctuation theory. If we consider the case with a small 
number density difference 5n/n = 0.1, we hnd that on 


Figure 3: (a) Tc as a function of Sn/n, with fixed 1/kFa = 
—0.5 and three different SOC coupling A = 0.25, 0.5, Ivf- 
(b) Tc as a function of Sn/n, with fixed 1/kFa = 0 and three 
different SOC coupling X/vf = 0.25, 0.5, 1. 


the BCS side the Tc is just slightly lower than that bal¬ 
anced case. But in the BEC side, this small number 
difference makes the small local minimum in the near 
BEC side disappear. Therefore, we are left with a little 
more broad maximum peak around the unitary limit. In 
the near BEC regime, the imbalanced Fermi gases have 
higher Tc than the balanced case, in contrast to the ex¬ 
pectation that the spin polarization should suppress su¬ 
perfluidity. This phenomenon actually already happened 
for very small Sn in the Fermi gases without SOC. The ef¬ 
fect of SOC is to make this non-monotonic effects appear 
in a much larger parameter space than before. If we tune 
the number density difference to dn/n = 0.3, then the 
whole Tc curve is evidently lower than the balanced case. 
We also note that the location of the maximum is shifted 
to the BEC side which is same as the Fermi gases with¬ 
out SOC. For the 5n/n = 0.3 case, our Tc curve stops 
at shallow BCS region because the un-physical behav¬ 
ior appears in a larger region due to the larger number 
difference. In panel (b) of FiglU we show the Tc curve 
with hxed number density difference 6n/n = O.I. In the 
case without SOC, there appear the un-physical reen¬ 
trant effects around the unitary and the BCS side as in 
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balance. 

In FigUl the as a function of \/vp is plotted, with 
different number density imbalance 5n/n = 0, 0.1, 0.4, 
respectively. The results of shallow BCS with 1/kpa = 
—0.5 and unitary limit are presented in panel (a) and 
(b) respectively. Since the SOC increases the pairing 
strength, the Tc v.s. X/vp plot can also be treated as 
another type of BCS-BEC crossover, with the large SOC 
corresponding to the deep BEC limit. When i5n = 0, 
we reproduce the results in [2^. One can see that Tc 
almost did not increase with the increasing of the SOC 
coupling. This is because we choose to study with shal¬ 
low BCS or unitary limit in order to avoid un-physical 
behavior due to the imbalance. Close to the unitary limit, 
Tc is already very high, thus obscure the effects of SOC. 
One can see that the Tc v.s. X/vp curve generally as¬ 
sume the similar shape as the Tc v.s. 1/kpa plots. The 
Tc reaches the highest value near the unitary limit then 
there is a small dip as one further increase the pairing 
strength. We should mention that there are some un¬ 
physical re-entrant behaviors appear in the small SOC 
area with large imbalance. Other than these small area, 
other part of the plot represents the stable polarized su¬ 
perfluid phase. 

From the numerical results showing above, we can con¬ 
clude that the SOC weaken the effect of the population 
imbalance. 


Figure 4: (a) Tc as a function of X/vf, with fixed l/kpa = 
—0.5 and three different number density difference 5n/n = 
0, 0.1, 0.4. (b) Similar as panel (a) but with interaction 

l/kpa — 0. 

[s^. However, when we turn on the SOC term, this un¬ 
physical behavior disappears. One can see that the SOC 
greatly stabilizes the polarized superfluid phase. 

In FigI3J we show the TcS as a function of num¬ 
ber density difference dn/n, for different SOC couplings 
X/vp = 0.25, 0.5, 1. The results of shallow BCS with 
1/kpa = —0.5 and unitary limit are presented in panel 
(a) and (b) respectively. In both cases, the increasing of 
5n will rapidly suppress Tc as one expected. As we dis¬ 
cussed before, the effects of SOC is to enhance the pair¬ 
ing between fermions thus make the system more BEC 
like. Here we only present the results with large enough 
SOC such that there is no un-physical behavior in the 
Tc curve. In the shallow BCS case, one can see that the 
increasing of A indeed push the whole Tc curve upwards. 
On the other hand, we note that in the unitary limit, the 
SOC term actually make Tc lower for small imbalance 
and increase Tc for large imbalance. The overall effect 
of SOC in unitary is making the Tc curve more flat than 
before. Similar behaviors are also observed in the BEC 
side. Therefore, for the unitary limit and the BEC side, 
the enhancement of pairing due to SOC roughly cancels 
out the suppression of pairing due to the population im- 


V. CONCLUSION 

In this paper, we have studied the BCS-BEC crossover 
of three dimensional spin polarized Fermi gases with 
Rashba SOC. To capture the fluctuations due to the 
non-condensed pairs at finite temperature, we introduce 
a single particle self-energy dressed by the T-matrix or 
pair propagator approximated by GqG ladder diagrams. 
Based on this GqG pair fluctuation theory, it is relatively 
straightforward to determine the pseudogap due to the 
stronger than BCS pairing. In the same time, the van¬ 
ishing of the order parameter determines the Tc, which 
is beyond the simple mean field approximation. 

With both SOC and spin polarization, there is a very 
large possible parameter space which may support a va¬ 
riety of exotic phases. In this paper, we focus on the 
simplest possible phase, the Sarma state or the spin po¬ 
larized superfluid state, which is a direct generalization of 
BCS superfluid state. First, we have shown that the be¬ 
haviors of the total gap and pseudogap with SOC are very 
similar to that without spin polarization, which suggests 
that the SOC increases the pairing effects and stabilize 
the spin polarized superfluid state. 

We then present the detailed results of various Tc 
curves which can also be thought as preliminary phase di¬ 
agrams. One can see that without SOC, the un-physical 
re-entrant behavior appears in the shallow BCS side, 
which signals the phase separation between BCS super¬ 
fluid state and spin polarized normal state. With certain 
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amount of SOC coupling, one can check that spin po¬ 
larized superfluid state are stabilized against the phase 
separation. From the Tc v.s. 5n curve, one can see that 
the Tc curves are flattened by increasing the SOC cou¬ 
pling, which implies that the effects of spin polarization 
and SOC largely canceled out. In summary, our work 
suggests that the stable region of spin polarized super¬ 
fluid state are greatly enlarged by the introducing the 
Rashba SOC term. 

The authors are supported by NSFC under grant No. 
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Appendix A: detailed expressions for constant ao 
and Bi 

In the appendix, we show the detailed expressions for 
the spectral weight factor ag and inverse pair masses Bi 
as follows. 
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where i = x,y and 


yl = S7(a^^ + a7^k£^k+Pk), B = \^(kl + kl) 

c = Pk^^k + 7Pk4 + a7CkPL D = S(5p^[7Pk + a7(A2 + ^l) + a^i,E^] 

J = — - L = AB-D-Crn, 

m 77k 

Ck-q = ~^k-q + S77k-q (^4) 
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